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ABSTRACT 

We compare the radio and soft X-ray brightness as a function of position within the 
young supernova remnant Cassiopeia A. A moderately strong correlation (r = 0.7) was 
found between the X-ray emission (corrected for interstellar absorption) and radio emis- 
sion, showing that the thermal and relativistic plasmas occupy the same volumes and are 
regulated by common underlying parameters. The logarithmic slope of the relationship, 
ln(Sx-ray) = 1-2 x l^SVadio) +ln(fc) implies that the variations in brightness are primarily 
due to path length differences. The X-ray and radio emissivities are both high in the 
same general locations, but their more detailed relationship is poorly constrained and 
probably shows significant scatter. 

The strongest radio and X-ray absorption is found at the western boundary of Cas A. 
Based on the properties of Cas A and the absorbing molecular cloud, we argue that they 
are physically interacting. 

We also compare ASCA derived column densities with A21 cm H I and A18 cm OH 
optical depths in the direction of Cas A, in order to provide an independent estimate 
of ISM properties. We derive an average value for the H I spin temperature of ~ 40°K 
and measure the ratio OH/H 2 , which is nominally larger than previous estimates. 

Subject headings: ISM: clouds, atoms, molecules, supernova remnants, individual: Cas A 
— radio lines: ISM — X-rays: ISM — radio continuum: ISM 
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1. Introduction 

The basic hydrodynamical structure of idealized 
young supernova remnants (SNRs) seems well-under- 
stood from a theoretical standpoint (Gull 1973a, Che- 
valier 1982). We expect to find an outer shock, a 
contact discontinuity between the shocked circumstel- 
lar medium and the ejecta, and a reverse shock mov- 
ing into and decelerating the ejected material. Each 
of these should give rise to radio and X-ray radia- 
tion, with different emissivities depending on the local 
physical processes. 

Obscrvationally, the situation is far from this ideal. 
In the best studied young SNR, Cas A, none of these 
structures can be clearly identified (Anderson & Rud- 
nick 1995, hereafter A&R). There are also questions 
about the nature of the outer shock, where the ex- 
pected tangential magnetic fields arc not seen (e.g., 
Kepler - Dickel et al. 1989, and Tycho - Dickel et 
al. 1991). Inhomogeneities on a variety of scales also 
complicate the observational as well as the theoreti- 
cal pictures (e.g. Borkowski et al. 1992, Cliffe & Jones 
1994, Jun & Norman 1994). 

We need to clarify the nature of actual SNR struc- 
tures both to understand the hydrodynamics and also 
to begin addressing important physical issues such as 
magnetic field amplification and relativistic particle 
acceleration. Although reasonable theoretical mecha- 
nisms exist for these processes (Gull 1973b, Reynolds 
& Ellison 1992), the observational signatures are far 
from clear (e.g. Anderson et al. 1994). One fruitful 
approach to addressing such questions may be a care- 
ful examination of the relationship between the X-ray 
and radio emissivities within a remnant, because of 
the different physical processes involved. 

The bulk of the X-ray emission at low energies re- 
sults from thermal line emission (Becker et al. 1979, 
Holt et al. 1994, hereafter HGTN), depending pri- 
marily on the temperature (« 3 keV) and density 
(« 10 cm~ 3 ) of the plasma carrying most of the mass 
and momentum (Fabian et al. 1980, hereafter F80). 
On the other hand, the radio emission is synchro- 
tron radiation from relativistic electrons (e « 0.05 — 8 
GeV) in magnetic fields of w 100 [iG (Cowsik & 
Sarkar 1980). In one remnant, SN1006, X-ray syn- 
chrotron radiation is probably present at keV energies 
(Koyama et al., 1995), although this is an exceptional 
case. 

Very little quantitative work has been done on 
the comparison of X-ray and radio emissivities in 



SNRs. The canonical wisdom is that the two are 
well-correlated on large scales, but show little correla- 
tion at smaller spatial scales (F80, Matsui et al 1984, 
hereafter MLDG). MLDG studied these relations in 
Kepler, where they divided the remnant into twelve 
sectors and found a moderate correlation of the form 
In (fix-ray) ~ (1-1 to 2.5) xln(S ra di )+ln(scale factor). 

We chose Cas A for study because of the availabil- 
ity of both high quality radio and X-ray data. In the 
cm wavelength range, Cas A is the brightest object in 
the sky outside of the solar system. At an age of 300 
years, it is believed to be in a pre-Sedov phase, and is 
situated 3.4t°; 3 kpc away (Reed et al. 1995), at the 
far edge of the Perseus arm. The column densities of 
hydrogen between here and Cas A (A H ~ 10 22 cm~ 2 ) 
are such that the optical depths of 1-2 keV X-rays 
are of order unity. Therefore, column densities inho- 
mogeneously distributed across the remnant, such as 
due to structures local to Cas A and the Perseus arm, 
can play a large role in determining the apparent soft 
X-ray morphology. For these same reasons, Cas A is 
an excellent choice as a background source for radio 
and X-ray interstellar medium (ISM) studies. 

The ISM has been well-studied in A21 cm absorp- 
tion (Mcbold and Hills 1975, hereafter MH75; Bieging 
et al. 1991, hereafter BGW; and Schwarz et al. 1996, 
hereafter SGK). Although BGW's VLA study was 
the highest resolution and most detailed, it only cov- 
ered the velocity range of the Perseus arm. SGK's 
Westerbork study, though only at a resolution of 30", 
covered both the Orion and Perseus spiral arms. In 
addition there have been numerous molecular absorp- 
tion studies using Cas A — in H 2 CO (Goss et al. 
1984), in CO (Troland et al. 1985, hereafter TCH and 
Wilson et al. 1993, hereafter WMMPO), NH 3 (Batrla 
et al. 1984 and Gaume et al. 1994) and OH (Bieg- 
ing & Crutcher 1986, hereafter BC). The absorption 
patterns in the various molecules are similar to each 
other, but very different than that of the H I. How- 
ever, C II seems to be correlated with the H I instead 
of the molecules (Anantharamaiah et al. 1994). 

Rasmussen (1996, hereafter R96) studied the spa- 
tial dependence of X-ray model parameters using the 
ASCA satellite, resulting in a total column density 
(Ah) map. In this paper, we compare the radio ab- 
sorption data of BC and SGK to the Ah map of R96. 
From this we measure the scaling relation between 
column density and equivalent line widths from the 
radio absorption measurements of SGK and BC. This 
allows measurements of the average H I spin temper- 
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ature and the Aoh/Ah 2 abundance of the ISM to be 
calculated. 

2. Analysis 

2.1. Comparing Cas A's X-ray and Radio 
Surface Brightnesses 

In order to compare the HRI X-ray (S'hri) and 
radio (SVadio) surface brightnesses, we must first cor- 
rect for the soft X-ray optical depths, which are of 
order unity. The higher energy ASCA data do not 
suffer from this problem, but a full analysis of those 
data at sufficient resolution has not yet been done. 
In order to estimate the X-ray optical depths, we 
made use of BC's A18 cm OH and SGK's A21 cm 
H I absorption data (see Fig. A & B and Fig. 0) 
to calculate the total column density (iVn, Fig- 
and Fig. ||) of hydrogen along each line of sight to 
the remnant. Once an Ah map has been calculated, 
a corrected image of unabsorbed soft X-ray emission 
(Sx ctd , Fig. [T]e) can be derived. Next, we can compare 
the actual thermally produced X-ray surface bright- 
ness (Sx ctd ) with Cas A's well-known radio synchro- 
tron morphology (SVadio, Fig. |l|F). We first describe a 
direct method for carrying out these corrections, and 
then the alternative scheme we found it necessary to 
use. 

We start by writing the X-ray optical depth as 

T X EE Of* X Nu (1) 

where c™ 1 1S ^ ne cross-section per hydrogen atom av- 
eraged over the HRI bandpass; Ah is the total num- 
ber of hydrogen atoms along the line of sight. We can 
then calculate 

A H = N m + 2 x A H2 (2) 
For the atomic hydrogen, 

A HI = (1.83 x 10 18 EW m ) (^f) cm" 2 (3) 

where T sp - ln is the spin temperature of the atomic hy- 
drogen, and the equivalent width in H I is calculated 
by integrating the observed optical depths t v over all 
velocities 

EW m = Jr v dv (4) 

In the case of the molecular hydrogen, for which direct 
measurements are not available, we can use the OH 




Fig. 1. — 30" resolution images of Cas A. The 
histogram-equalization method of scaling was used to 
enhance the images. The quantities represented are: 
(A) the equivalent width of the A21 cm line; (B) the 
equivalent width of the A18 cm (OH) line; (C) the 
total column density as derived from images A and 
B; (D) the logarithm of the ROSAT HRI image; (E) 
the logarithmic HRI image corrected for absorption; 
(F) a logarithmic A20 cm VLA continuum map. His- 
tograms of the quantities shown in images A, B and C 
are shown in figure ||; image D ranges from -2.6 - -0.8 
ln(cts s _1 beam -1 ); a plot of image E versus image F 
is shown in figure m 




H I Equivalent Width ((In/,)-') OH Equivalent Width ((km/.)"') Calunnr. Dan.it, (10= cm"') 



Fig. 2. — Histograms of A21 cm equivalent width, A18 
cm equivalent width and the derived column density 
{N u = D(EW m ) + E(EW n)+F). 
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absorption data, and write 



A^OH 



2.2 x 10 14 EWon) 



T 

ox 

°K 



(5) 



with equivalent width defined as above. 

This approach (of finding the X-ray absorbing col- 
umn densities a priori from radio line data) is the 
straight forward one, but it has several problems be- 
cause of uncertainties in the scaling constants that 
must be used. The effective X-ray cross-section has 
been modeled by Morrison and McCammon (1983, 
hereafter M&M), who showed that the monoener- 
getic cross section depends strongly on photon energy 
and changes discontinuously at quantum thresholds. 
Therefore, ct™ 1 is highly dependent on Cas A's soft 
X-ray spectrum, and thus has significant uncertain- 
ties. In addition, it is difficult to know what value 
to use for the H I spin temperature, as was demon- 
strated by MH75. Another major uncertainty comes 
from the OH abundance value. These issues will be 
discussed more fully in section 
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We therefore adopted a different approach, simul- 
taneously determining the X-ray/radio relationship 
and the absorption corrections by minimizing the 
quantity 



X 



Ms* 



radii 



ln(fc)] 5 



+ {r/ i _ 1 X (Tradio) 5 



(G) 



with respect to the parameters i], A, B, and [C — 
ln(fc)], where the logarithmic corrected X-ray image 
is given by: 



ln (Sx ctd ) = 



ln(5 H Ri) + t x 
ln(5 HR i)+Ax EW m 
+B x EWqh + C . 



(7) 



The parameter rj measures the logarithmic scaling 
between the X-ray and radio emission. The parameter 
[C — ln(fe)] contains both information on the normal- 
ization of emissivities (equation^, parameter \n(k)), 
as well as allowing for absorption that is not well- 
modelled by the H I and OH optical depths (equa- 
tion fj], parameter C). This could be due, e.g., to H i 
saturation, variations in T sp ; n , which are not reflected 
in X-ray absorption, and the likely existence of Orion 
Arm molecular gas outside of the BC measurements. 

Because the errors in the data are a function of 
the radio and X-ray flux, it is important that the er- 
ror images are included in the determination of \ 2 
(equation ||) . 



As a simplifying procedure, we performed the fits 
iteratively, holding the errors fixed as propagated 
from the previous best fit parameters (A;_i, 
and ?7 i _ 1 ): 



CT ln(S H 



hi) + (^-! x a Ew m y 

(-Bj-i x &ew h) 2 ■ 



(8) 



Prior to performing the x 2 analysis, the maps were 
prepared as follows. In order to match the resolu- 
tion of the H I data, we smoothed the OH line and 
continuum maps to a resolution of 30", before cal- 
culating the optical depths and equivalent widths. 
We also smoothed the epoch 1990 (3.4 ks live-time) 
ROSAT HRI image (obtained from the HEASARC 
public archive) and the epoch 1987 A20 cm VLA con- 
tinuum map described by A&R, to 30" resolution. 
Errors in each map were calculated using standard 
methods. The resulting maps contain 51 independent 
beams, and are shown in figure [j]. Histograms of the 
radio absorption values are shown in figure |[ 

2 

The minimum value of jr was 47, showing that 
there is a large amount of scatter still unaccounted for 
in the X-ray/radio relation. The number of degrees of 
freedom, N, coincidentally, also equals 47. In order to 
calculate 70% and 99% confidence limits (5xtq,5x99) 
for parameter X , we determined the value of 



X 2 (X±5x70. 



X99) 



Xn 



(1.097,1.542) (9) 



where Xmin ^ s the global minimum value of \ 2 - Fol- 
lowing Bevington & Robinson (1992), while parame- 
ter X is being varied, all other parameters are allowed 
to float to minimize % 2 for that value of X. Figure || 
shows contour plots of the confidence levels, as a func- 
tion of parameters A and B, i] and [C — ln(fc)]. 

Figure^ shows the correlation between the thermal 
X-ray emission ln(5x ctd ) an d radio synchrotron radia- 
tion (ln(Sradio))- The best fit slope is ry = 1.21 ±0.41, 
(70% confidence limit), with the error probably dom- 
inated by real scatter in the X-ray/radio relation. In 
section ||, we will discuss the competing effects of 
emissivity and path length that can contribute to such 
differences. 

2.2. Interstellar Medium Parameters 

In this section, we introduce a technique for de- 
termining physical parameters of the ISM. Although 
our results have large uncertainties, the method is 
useful in confirming previous estimates, and could be 
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Fig. 3. — 70% and 99% confidence levels derived from 
the goodness of fit x 2 obtained by fitting a power- 
law relation between the radio and corrected X-ray 
morphological distributions. They are shown here as 
a function of the parameters A, B, [C — ln(fc)] and r) 
as defined in the text. The x in each plot represents 
the minimum \ 2 ■ 



I + 




3.0 - - 



2.7 3.0 3.3 3.6 3.9 

ln(S„„„) (Jy/Deam) 

Fig. 4. — Plot of the logarithmic soft X-ray emis- 
sion corrected for absorption using the best fit val- 
ues for parameters A, B, and [C — ln(fc)] as defined 
in equation versus the logarithmic A20 cm emis- 
sion. The line represents our the best fit model: 
ln(Sx ctd ) - ln(fc) = r?ln(S ra dio). 



extended to higher accuracy. Using the above pro- 
cedure, we simultaneously determined the scaling pa- 
rameters between H I and OH optical depths and soft 
X-ray absorption due to the line-of-sight ISM column 
density. The best fit values are A = 0.073 ± 0.039 
(km/s)-\ B = 1.6 ± 1.0 (km/s)- 1 and [C - ln(fc)] = 
1.6 ±2.2 (Fig. |). 

The X-ray absorption and the total column den- 
sity can also be measured more directly through X- 
ray spectral fitting. With spectral fitting, the effec- 
tive X-ray cross-section, which changes strongly as 
a function of photon energy (M&M), can be explic- 
itly included. R96 assumed Cas A's X-ray spectrum 
to consist of power-law emission plus many Gaussian 
spectral emission lines (see HGTN), absorbed by the 
ISM assuming M&M's effective cross-sections. This 
resulted in a 32 x 32 map (19 "/pixel) for each param- 
eter fit, including the column density. A factor of 2 
oversampling was used, so the "beam size" is approxi- 
mately 40"; this is also a function of ASCA's complex 
point spread function (Tanaka et al., 1994). However, 
since these fits did not cover the whole image of Cas 
A, and because they rely on a "super-resolution" of 
the data, we chose not to use them for the original 
correction of the X-ray maps. We will, however, use 
the ASCA data to derive parameters of the interstel- 
lar medium. 

In this section, we compare our optical depth im- 
age to R96's ASCA iVu image to derive a second set of 
scaling parameters for EWm and EWon ■ Combining 
this with the previously determined scaling parame- 
ters, and using equations || and ^, we can then deter- 
mine both the effective ROSAT HRI cross section as 
well as some important ISM parameters. 

To perform this analysis, we first aligned the ASCA 
field with our radio images by maximizing the cross- 
correlation of R96's -/Vh asca image with our optical 
depth image (A x EWm +B x EWqb)- This yielded a 
total of 25 independent 30" beams with which to com- 
pare the X-ray and radio absorption. With -/Vh asca 
in real physical units (10 22 cm -2 ) and the equivalent 
widths in km/s, we then minimized the quantity: 

X 2 (D, E,F)^J2 ^Hasca ~ ^H Radlo ) 2 (10) 

where D, E, and F are defined by 

^H Radlo = D (EW m ) + E (EWoh) + F (11) 

The best fit values for parameters D,E and F are: 
(7.1±3.8) x 10~ 3 , 0.213±0.141, and 0.69±0.19 respec- 
tively (70% confidence); the x 2 dependence on D,E 
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and F is shown in Figure ||. The confidence limits are 
defined as: 



X 2 (X ± 6 X70 , xg9 ) = (1.134, 1.831) x x 



2 

min 



(12) 



as appropriate for 22 degrees of freedom. The use 
of Xmin i s necessary because R96 did not provide an 
error image. 

At this point, we have two sets of scaling param- 
eters for EWni and EWon- One set depends on the 
HRI cross section; the other does not. We can there- 
fore use equations |?] and |ll| to find that 

A B 

°™ = D = E < 13 > 

This yields % = 10.3 ± 7.7 and § = 7.6 ± 6.8 with 
a weighted average of cr" RI = (8.8 ± 5.1); all are in 
units of 10 _22 cm 2 . This value of cr™ is several times 
larger than the rough estimate we made by averag- 
ing the M&M cross-section values over the bandpass, 
weighted by the HRI effective area, and assuming a 
flat spectrum for Cas A. If this discrepancy is real, 
then it indicates that on average, Cas A's source spec- 
trum is probably decreasing with increasing energy in 
the HRI bandpass. 

From our measured parameters D and E, and as- 
suming that we have isolated the contributions due 
to H I and OH in our x 2 analysis, we can then use 
equations || and || and the values of D and E to cal- 
culate T sp i n and Non/N-a 2 . However, the presence of 
a significant non-zero value for F shows that we have 
not successfully modelled all of the X-ray absorption, 
and we will have to take this into account. 

First, we consider the molecular hydrogen compo- 
nent. Figure ^ shows that there is very little corre- 
lation between parameter E and either parameter D 






Fig. 5. — The x 2 between the ASCA and radio ab- 
sorption line derived column densities as a function of 
the scaling parameters D, E and F, which are defined 
such that Nh = D (EW m )+E(EWon)+F. At each 
point in the plots, the parameter not shown was set 
to the value which would minimize x 2 there. Contour 
levels represent confidence levels of 70% and 99%. 



or F. This implies that our best fit parameter E is 
robust to uncertain contributions to the absorption 
from other ISM components. We therefore find the 
fraction of OH in the Perseus spiral arm as: 

t: = (4.1 ± 2.7)xl0-=(^). (14, 

This derived ratio is nominally higher than the (2.9 ± 
2.7) x 10" 7 quoted by WMMPO. 

The contribution of the atomic component (param- 
eter D) is more problematic, because of its correlation 
with parameter F (see Fig. ||). Formally, we can cal- 
culate an average spin temperature for the H I as 

T spin = D/(1.83 x 10~ 4 ) = (39 ± 21) °K . (15) 

This temperature is the same as the spin temperature 
derived by Payne et al. (1994) towards Cas A and in 
the range of those found by Kalberla, et al. (1985), 
toward 3C147. The major contribution to the error 
in D actually comes from its correlation with F. 

There are three likely sources for a non-zero F: 
Orion arm molecular gas (not measured by BC), sat- 
urated H I and hot intercloud H I. 

TCH observed Orion arm 13 CO column densities 
to be approximately 20% of the Perseus arm column 
density. This implies that only about 6 x 10 20 cm 2 
of -/Vh asca (» 10% of parameter F) can be accounted 
for by molecular gas outside of BC's velocity range. 

BGW named the highly saturated H I absorption 
feature near a velocity of -48 km/s "the curtain", be- 
cause of its spatial uniformity and high optical depth 
(mostly above 5). Here SGK's H I equivalent width 
measurements underestimate the actual absorption. 
If SGK accounted for about half of the H I in "the 
curtain" , the unmeasured portion should account for 
about 1-2 x 10 21 cm 2 of Na ASCA or « 15% of parameter 
F)- 

MH75 used the Effclsbcrg 100m telescope to ob- 
serve the A21 cm line in both emission and absorption 
toward Cas A. Their observations are consistent with 
a two-temperature H I model; where the hot compo- 
nent has spin temperatures in excess of about 3500 
°K, which would be observed only minimally by SGK 
(equation ||). MH75's estimated cloud/intercloud 
mass ratio is about 1:2, which could account for all of 
our measured parameter F. 

These factors complicate our determination and in- 
terpretation of T sp i n , as described above, so it should 
be considered simply as a characteristic value for the 
cooler gas. 
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3. Discussion 

3.1. The X-ray/Radio Emissivity Relation 

The long-term objective of comparing the X-ray 
and radio emission in Cas A is to identify the state 
and structure of the thermal and relativistic plasmas 
and the physical connections between them. The X- 
ray emission spectrum from Cas A is itself believed 
to come from two different plasmas (Janscn ct al. 
1988). The low temperature plasma (kT < 1 keV) 
is most likely reverse-shocked ejecta; models of the 
emission in the ROSAT band show that the emissivity 
is a mixture of thermal bremsstrahlung and line emis- 
sion. The higher temperature plasma (kT sa 3 keV) is 
most likely circumstellar matter shocked by the blast 
wave. The reverse shock in Cas A is believed to cur- 
rently dominate the soft X-ray emission (F80, Jansen 
et al.); the factors that influence the relative lumi- 
nosity of these two shocks are discussed in detail by 
Masai (1994). We have shown here the importance of 
absorption for the low energy X-rays; analyses based 
on hardness ratios (e.g., HGTN) must first correct for 
this effect. 

In the radio band, most of the emission comes from 
the bright ring, which is identified with either the re- 
verse shock, the contact discontinuity, or both (Gull 
1973 a, b). There is also a lower surface brightness 
radio plateau beyond the ring and a wealth of struc- 
tures on smaller scales, including knots, bow shocks, 
filaments, etc, with lifetimes of order 30 years. 

To interpret the observed radio/X-ray correlation, 
we must now distinguish between surface brightness 
and emissivity. If the emissivities in a remnant do not 
vary spatially, then all variations in surface brightness 
must be due to variations in path length (or filling 
factor), and the X-ray and radio brightnesses should 
be proportional to each other (logarithmic slope of 
1). Some of the brightness variations in Cas A are 
clearly due to path length differences, such as the 
bright ring itself. However, even within the bright 
ring, the local radio brightness may largely be due to 
the path length through the emitting material at that 
position, rather than large variations in the intrin- 
sic emissivity. The observations of Reed et al. (1995) 
show that the apparently complete Cas A shell is very 
non-uniform in optical line emission. The analysis of 
A&R show the presence of major dynamical asymme- 
tries which are probably coupled to spatial variations. 
The large velocity gradients in the X-ray data from 
ASCA (HGTN) also demand large-scale non-spherical 



structures in the X-ray emitting material. 

Given the major role played by path length vari- 
ations through the emitting material, our observed 
logarithmic slope of 1.2 for the X-ray/radio surface 
brightness relationship can be understood only as a 
lower limit to the slope relating the actual emissiv- 
ities. If the variations in emissivity are much less 
than the variations in path length, the emissivity 
slope could be much higher. We are therefore unable 
to comment on how the various processes leading to 
soft X-ray and radio emission vary with one another. 
The same problem might easily affect the results of 
MLDG, whose X-ray/radio logarithmic slope could be 
as low as 1.1. At present, it is not clear how to isolate 
true variations in emissivity. 

The modelling of the radio/X-ray emissivity rela- 
tionship in an inhomogeneous rapidly-evolving rem- 
nant such as Cas A is also quite uncertain. It is im- 
portant to avoid the simple scaling relations based on 
X-ray thermal bremsstrahlung, such as discussed by 
MLDG, because such relations ignore the dependence 
of line emission on factors other than density. In ad- 
dition, we know from recent numerical simulations 
that the magnetic field amplification (Jun & Norman 
1994) and relativistic particle acceleration (Jones et 
al. 1994) reflect the history of the plasma and cannot 
be simply described by current state parameters such 
as density and temperature. Jun (1995) has mod- 
eled the radio synchrotron radiation and the thermal 
bremsstrahlung component of X-ray emission in his 
3D MHD simulations of young supernova remnants. 
Although he finds the same major features in the X- 
ray and radio emission, e.g., the clumpy bright ring, 
there is only a weak correlation between the two (Jun, 
1996). This is due to the strong dependence of only 
the radio emission on the local magnetic field. There- 
fore, in order to effectively use such observations as 
presented here and in MLDG, we need both more 
sophisticated time-dependent analyses of the X-ray 
radiative transfer and an understanding of the rela- 
tionship on various scales between the magnetic field 
and other hydrodynamical parameters. 

3.2. The Western Molecular Cloud 

The western edge of Cas A is unusual in a num- 
ber of ways, suggesting an interaction between the 
expanding SNR and a local molecular cloud. Con- 
sidering the properties of Cas A itself, we first note 
that the brightest radio and X-ray emission is found 
in this region. HGTN show that the western region 
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has significantly lower X-ray equivalent widths than 
the other major emission regions. They interpret this 
as due to relatively stronger emission from the outer 
shock, which would be expected if it were moving into 
an area of higher density. A&R found that the mo- 
tions of bright radio knots in the west showed extreme 
departures from the relatively uniform expansion seen 
in other regions. Many western knots are actually 
moving back toward the center of expansion - these 
must be due to an external interaction. This is also 
a region of steep radio spectral indices, implying that 
the conditions for relativistic particle acceleration are 
different here (Anderson & Rudnick, 1996). 

On the opposite side of the remnant, there is a 
break in the shell, and groups of fast moving radio 
(A&R) and optical (Fesen et al., 1988) knots. This so- 
called "jet" could result from expansion into a lower 
density region opposite the molecular cloud. A similar 
situation is modelled by Tenorio-Tagle et al. (1985). 

Turning now to the external material, figure |l| (im- 
age C) is the derived Nn map, which shows strong 
absorption on the western side of Cas A. This dense 
cloud also shows up as the extreme value in the col- 
umn density histogram in Fig. ||. The cloud is at an 
LSR velocity of « -40 km/s, placing it in the Perseus 
arm (Batrla et al. 1984, Goss et al. 1984, TCH, BC 
and WMMPO). The OH column density maps of 
BC around -40 km/s trace out the regions of anoma- 
lous radio knot velocities and steep spectral indices 
discussed above. A trace of this cloud may also be 
visible in the H I measurements of BGW and SGK. 

Gaume et al. (1994) studied the NH 3 and CO ab- 
sorption towards the bright western region. For the 
-39 km/s cloud, they determined the density to be 
«h 2 ~ 1000 cm -3 characteristic of dark dust clouds, 
but found a higher than average kinetic temperature 
(m 18 °K) and line width (3.5 km/s). They suggested 
that the high temperature could be due to either an 
increased level of cosmic rays, or by cloud-cloud colli- 
sions. On the basis of both the unusual properties of 
Cas A in the west and the unusual cloud conditions 
there, it thus seems quite likely that an interaction is 
currently in progress. 

Wilson and Mauersberger (1994) pointed out that 
the circular shape of Cas A argues for a lack of ex- 
ternal interaction. However, the bright radio ring 
illuminates material that has just recently been de- 
celerated, because it is found in the same area as the 
optical fast moving knots, which are travelling several 
times faster (A&R). R. Dohm-Palmer & T.W. Jones 



(1996) have performed 2D numerical simulations of 
an SNR expanding into a sharp ISM gradient. They 
find that at the time when the reverse shock on the 
high density side has reversed the velocity of some 
compact features, the overall ring deviates from be- 
ing circular only by « 10%. This is the same degree of 
non-circularity observed for Cas A's bright ring, and 
so removes this objection to an external interaction. 

4. Conclusion 

In this paper we have presented a technique to cor- 
rect for spatially inhomogeneous absorption of soft X- 
rays in Cas A using radio absorption data. We find 
a good correlation between the soft X-ray and radio 
synchrotron emission from Cas A, but with signifi- 
cant scatter. The correlation is probably dominated 
by variations in path length, implying that the X-ray 
and radio emissions both occupy the same volumes. 
However, we have no evidence for a more detailed 
relation between their emissivities. A quantitative 
interpretation of these results requires more sophis- 
ticated modeling of both the X-ray radiative transfer 
and the relativistic plasma evolution in young SNRs. 

Future X-ray/radio comparisons of Cas A should 
concentrate in at least the following two directions: 
studies at higher spatial resolution with deeper HRI 
measurements and comparisons with ASCA's spa- 
tially resolved spectroscopy. With a deeper ROSAT 
HRI observation, it may be possible to separate emis- 
sivity from path length variations. In addition, X- 
ray proper motions could be measured and compared 
with the radio proper motions. Since HGTN's paper, 
Cas A has been used as a calibrator for ASCA, so 
more ASCA data have been obtained and ASCA's re- 
sponse functions have been refined. This will enable 
image reconstruction techniques to be applied to nar- 
rower bandpass images and better quality spatially 
resolved spectral fitting. Studying the correlations 
between the radio and X-ray morphologies as a func- 
tion of X-ray energy will allow the different emission 
mechanisms and temperature and metalicity struc- 
tures to be distinguished. 

We have shown that Cas A is likely to be interact- 
ing with a dense cloud in the west. This has affected 
both the properties of the remnant and the cloud. 
Such interactions may play an important role both in 
SNR dynamics, and in the transfer of energy into the 
ISM. 

We have also demonstrated a new technique for 
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probing the ISM. By comparing X-ray and radio spec- 
troscopic absorption measurements, the H I spin tem- 
perature and molecular abundances ratios were mea- 
sured. Future studies of other radio and X-ray bright 
extended objects can significantly enhance our under- 
standing of the ISM, by comparing spatially resolved 
column densities from either the ROSAT PSPC or 
ASCA with radio and far IR atomic and molecular 
line data. 
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